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We analyze the conventional perturbative treatment of sphaleron-induced baryon number washout relevant 
for electroweak baryogenesis and show that it is not gauge-independent due to the failure of consistently imple- 
menting the Nielsen identities order-by-order in perturbation theory. We provide a gauge-independent criterion 
for baryon number preservation in place of the conventional (gauge-dependent) criterion needed for success- 
ful electroweak baryogenesis. We also review the arguments leading to the preservation criterion and analyze 
several sources of theoretical uncertainties in obtaining a numerical bound. In various beyond the standard 
model scenarios, a realistic perturbative treatment will likely require knowledge of the complete two-loop finite 
temperature effective potential and the one-loop sphaleron rate. 



(N 



Oh- 
I 

Oh- 

CD 



> 

in 



X 



I. INTRODUCTION 

Explaining the cosmic baryon asymmetry remains an open 
problem at the interface of cosmology with particle and nu- 
clear physics. Among the several proposed scenarios, elec- 
troweak baryogenesis (EWB) is particularly interesting as it 
is conducive to experimental tests. The mechanism requires a 
strong first order electroweak phase transition (EWPT) that 
proceeds via bubble nucleation at temperatures ~100 GeV 
and sufficient CP violation. CP-violating interactions at the 
bubble walls induce a net density of left handed fermions 
(n^) that biases electroweak sphalerons into generation of 
baryon number density. As the transition proceeds, baryon 
number diffuses into the interiors of the expanding bubbles 
where electroweak symmetry breaking slows the sphaleron 
transitions. The latter leads to baryon number erasure as the 
sphalerons try to restore chemical equilibrium. Thus, success- 
ful EWB requires that the sphaleron rate is sufficiently sup- 
pressed inside the bubbles to prevent this washout. 

In the standard model (SM), the effects of CP-violation as- 
sociated with the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix are too feeble to generate a sufficiently large ul to bias 
the sphalerons in the first place. Even if this were not the 
case, however, the finite temperature dynamics of the SM 
scalar sector do not allow for a first order EWPT as needed 
to prevent erasure of an initial baryon asymmetry. Conse- 
quently, successful EWB requires both additional sources of 
electroweak CP-violation and an augmented scalar sector. Ex- 
perimentally, searches for the permanent electric dipole mo- 
ments of the neutron, electron, and neutral atoms provide a 
probe of the requisite new CP violation, while searches for 
new scalar particles at the CERN Large Hadron Collider may 
reveal the ingredients needed for the strong first order EWPT 
that would preserve the resulting baryon asymmetry. 
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In this work we focus on the theoretical analysis of the 
washout of the baryon asymmetry during the phase transition. 
The most theoretically robust approach to the study of EWPT 
dynamics involves Monte Carlo lattice simulations. Several 
such studies have been caiTied out in the SM, yielding re- 
sults for the critical temperature, bubble nucleation rate, and 
sphaleron rate [1-7]. These quantities depend critically on the 
Higgs quartic self coupling, which also governs the value of 
the Higgs boson mass, wih- Consequently, there exists a tight 
connection between mu and baryon number washout. Lat- 
tice results indicate that a strong first order EWPT as needed 
for baryon number preservation requires mn ^ 70 GeV, well 
below the current LEP lower bound of 114.4 GeV[8]. Thus, 
EWB is only viable in the presence of an extended scalar sec- 
tor that is not subject to the current LEP limits. 

A number of extended scalar sector scenarios have been an- 
alyzed that may successfully remedy the SM shortcomings. 
The most widely considered are supersymmetric extensions 
(see e.g. [9-15] and references therein), though recent work 
has also included analysis of non-supersymmetric models (see 
e.g. [16-18] and references therein) . While there exist a 
handful of Monte Carlo studies of phase transition dynam- 
ics in the minimal supersymmetric standard model (MSSM) 
[19, 20], the majority of phase transition studies in extended 
scalar sector scenarios have relied on the use of perturbation 
theory. Moreover, non-perturbative calculations in the SM and 
MSSM have generally concentrated on the thermodynamic 
properties of a possible EWPT (such as the critical temper- 
ature) rather than directly computing the sphaleron rate that 
governs the degree of baryon number washout (for important 
exceptions, see [2, 7, 21]). 

This situation is not surprising, given the numerical cost 
in performing non-perturbative computations. At present, it 
is simply not feasible to survey the broad range of extended 
scalar sector models and to analyze in each the often large, 
multidimensional parameter space using Monte Carlo meth- 
ods. In order to identify the EWPT-viable parameter space re- 
gions and identify their phenomenological signatures for col- 
lider searches, the reliance on perturbation theory appears to 
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be unavoidable. Thus, it is desirable to arrive at the most 
theoretically sound perturbative treatment of baryon number 
erasure in order to identify candidate scenarios that are most 
promising for EWB and to determine their prospective exper- 
imental signatures. 

In what follows, we revisit the conventional, perturbative 
treatment of baryon number washout and observe that it suf- 
fers from a serious theoretical shortcoming, namely, depen- 
dence on the choice of gauge. Consequently, there exists rea- 
son to question conclusions drawn about the EWB -viability 
of scenarios analyzed according to the conventional perturba- 
tive treatments. We subsequently show that there exist means 
of arriving at a gauge-independent, perturbative treatment of 
baryon number washout and we analyze one such approach in 
detail. We argue that the use of this approach can provide a 
theoretically unambiguous indication of the conditions under 
which a given scenario may lead to baryon number preserva- 
tion and can point to the relevant, but more limited, regions 
of parameter space that should be explored with Monte Carlo 
computations. In the course of our discussion, we also re- 
visit the approximate criterion used to evaluate the efficacy 
of baryon number preservation — henceforth referred to as the 
"baryon number preservation criterion", or BNPC. We also 
argue that a realistic BNPC may be more stringent than the 
one widely used in recent years and that there exist important 
theoretical uncertainties that call for future study. As a practi- 
cal matter, we provide an approximate BNPC that differs from 
the standard formula used in the literature and that takes into 
account the foregoing issues. 



A. Brief Overview 

Before analyzing these considerations in detail, we provide 
a brief overview. If a first order EWPT occurs, it commences 
via the formation of broken electroweak symmetry bubbles at 
the nucleation temperature T/y. The latter typically falls just 
below the critical temperature Tc at which the broken and 
unbroken minima of the scalar effective potential become de- 
generate. Inside the broken symmetry phase, the depletion of 
baryon number density is expected to follow a first order 
rate law 



the "washout factor" 
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(1) 



where k{T) is the temperature-dependent rate constant[22] 
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governed by the temperature T, number of fermion gener- 
ations nf, and most importantly the sphaleron rate per unit 
volume Ts-ph/V that depends exponentially on the sphaleron 
energy relative to that of the electroweak vacuum, Ai?sph. Ad- 
ditional, non-exponential dependence on T is contained in the 
prefactor A{T). 

The extent to which an initial baryon asymmetry is erased 
over the duration of the transition, Atew, is characterized by 



n5(A^Ew) 
nB(0) 



(3) 



Assuming for illustrative purposes that k{T) is approximately 
constant during the transition, integration of ( 1 ) yields 



(4) 



where we have taken the critical temperature Tc as a char- 
acteristic temperature of the transition. The magnitude of an 
initial baryon asymmetry, determined by the strength of the 
CP-violating interactions and transport dynamics, will dictate 
a lower bound on S 



S>e- 
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below which the final baryon asymmetry would be smaller 
than the observed value. Baryon number preservation thus im- 
plies roughly a requirement on Ai^sph/Tc and thereby on the 
underlying particle physics that determines the critical tem- 
perature and the sphaleron energy at that temperature. 

In the conventional perturbative analysis, the sphaleron en- 
ergy Ai?sph is expressed in terms of (t)mm{Tc), the classical 
field that minimizes the finite temperature scalar effective po- 
tential at the critical temperature. Taking X Ri 10 (which we 
question below) and including the dependence on the At-Ew 
and other factors that enter the rate equation via (1) (see Ap- 
pendix A for details), one arrives at the following requirement: 



Tc 



> 1.0 (conventional criterion) . (6) 



As we discuss below, this criterion is open to question for 
three reasons: 

(1) The classical field 4>min{T) is inherently gauge- 
dependent at any temperature. 

(2) The method of computing Tc introduces additional 
gauge dependence that does not compensate the gauge- 
dependence of <f>mm{Tc)- 

(3) The choice of X ^ 10 is overly optimistic for most 
scenarios, and there exist additional uncertainties asso- 
ciated with computing the prefactor A{T) and integrat- 
ing (1). 

In what follows, we show that it is possible to obtain a 
gauge-independent perturbative estimate of the sphaleron rate 
and critical temperature and review both the theoretical uncer- 
tainties alluded to in point (3) as well as the choice of X. In 
brief, 

(1) The sphaleron rate T^ph/V is obtained by evaluat- 
ing the temperature-dependent effective action of the 
sphaleron. A gauge-independent, perturbative compu- 
tation of this quantity can be performed using a dimen- 
sionally reduced 3D effective theory where only the 
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gauge-independent ©(T^) terms are included in the ef- 
fective action. As a result, the theory contains a gauge- 
independent energy scale v{T) that is non-vanishing be- 
low a temperature To (different from Tc) and that char- 
acterizes the sphaleron energy. 

Although additional gauge-independent, one-loop ther- 
mal contributions to Fsph/T^ have been obtained numer- 
ically for the minimal SM [23-25], similar perturbative 
computations do not exist for its extensions. Conse- 
quently, the conventional practice has been to estimate 
these contributions by replacing the gauge-independent 
sphaleron scale v{T) with the gauge-dependent Higgs 
field 4>mm{T) that minimizes the full gauge-dependent 
effective potential Veff{4>, T). As a result, the standard 
practice leads to a gauge-dependent estimation of the 
sphaleron rate regardless of the temperature. Clearly, 
a gauge-independent calculation requires utilizing the 
gauge-independent scale v{T) rather than (t>mm{T)- 

(2) The critical temperature is obtained from the behavior 
of the finite-temperature effective potential 14ff((/),T). 
Bubble nucleation characteristic of a first order phase 
transition requires the existence of a barrier separating 
the electroweak symmetry-broken and -unbroken min- 
ima in Fcff (0, T). At one-loop order in the SM, this 
barrier is generated by non-analytic oc T4>^ terms in the 
potential. Since the early work of Dolan and Jackiw, 
however, it has been known that this non-analytic part 
of the potential is gauge-dependent [26], and the early 
papers argued that one should drop these terms from 
T4fi(0,r) when deriving Tc- Since doing so allows 
only for a second order phase transition, the common 
practice in EWB studies has been to retain the non- 
analytic terms in VJ;ff((/), T). This causes Tc to inherit 
the gauge-dependence. 

A gauge-independent value for Tc can be obtained in 
perturbation theory by either using a gauge-invariant 
source term in the generating functional or by 

working with a source term which not gauge- 
invariant, and consistently implementing the Nielsen's 
identity [27]. In this work, we follow the latter approach 
and defer an exploration of the former to subsequent 
study. Nielsen's identity, derived shortly after the work 
of Dolan and Jackiw, implies that the exact effective 
action Feff and potential V^ff are gauge-independent at 
their extremal points, but that mean fields (^min) derived 
from them remain gauge-dependent. Since Tc is de- 
fined in terms of minima of V^ff and the sphaleron rate 
is derived from the saddle point of Feff, both quantities 
are expected to be gauge-independent. However, naive 
truncation of the perturbative series introduces spurious 
gauge-dependence in Tc- As we show below, it is pos- 
sible to avoid this spurious gauge-dependence by deter- 
mining Tc from the extrema of 14ff((^, T) in a manner 
consistent with the ^.-expansion. 

(3) Provided that one obtains a gauge-independent com- 
putation of the sphaleron rate and critical temperature. 



there exist remaining uncertainties in using these quan- 
tities to derive a BNPC. For example, the choice ofX^ 
10 derives from early work by Shaposhnikov[28], who 
estimated the SM baryon asymmetry from the equation 
of motion for Chern-Simons number. The computation 
required an estimate of the quark-antiquark asymme- 
try (Sins — ^produced by CP-violation associated with the 
CKM matrix — that enters as a source term in the Chern- 
Simons equation of motion. While a range of possi- 
ble values for (5ms were considered in Ref. [28], the 
most generous possibility allows for a significant over- 
production of baryons in the SM and requires a large 
washout of the baryon asymmetry {X ^ 10). It was 
subsequently argued, however, that CP-violation from 
the CKM matrix would have been far too small to over- 
come even a washout factor of 0(1) (for a discussion, 
see e.g. Ref. [29]). 

Recent attention has focused on the production of CP- 
violating asymmetries in supersymmetric models, gen- 
erated by the scattering of superpartners from the bub- 
ble walls (see Ref. [30] and references therein). Con- 
straints on the supersymmetric CP-violating phases im- 
posed by EDM searches generally imply that these 
asymmetries are right on the edge, at best, of inducing 
the observed baryon asymmetry. Consequently, in these 
scenarios, a washout factor closer to 0(1) {or X ^ 1) 
is likely to be more realistic. In addition, the work of 
Ref. [7] suggests that at least in the Standard Model, 
the duration of the transition is likely to be closer to 
A^EW ~ 10^'^ti/, where in is the Hubble time, rather 
than A<EW ^ tn as assumed in arriving at the criterion 
in Eq. (6). 

There exist additional uncertainties associated with the 
frequency of the unstable sphaleron mode (known from 
the work of Ref. [23] but not generally included) and 
the evaluation of the sphaleron fluctuation determinant 
(known and generally considered) that should be taken 
into account when seeking to determine whether or not 
the underlying particle physics dynamics lead to suffi- 
cient preservation of an initial baryon asymmetry. The 
approximate BNPC that we derive below attempts to 
take these considerations fully into account. 

The primary implications of our analysis of these issues 
is to replace the gauge-dependent BNPC (6) with the gauge- 
independent one and to include the aforementioned theoretical 
uncertainties. In particular, in the LHS one should make the 
replacement 

ipmmiTc; v(Tc) 

TciO Tc ' 

where ^ is the gauge-fixing parameter and where the quanti- 
ties to the right of the arrow are ^-independent. On the RHS of 
the BNPC (6), one replaces unity by a range of bounds deter- 
mined by an appropriate choice of X and the other theoretical 
inputs. The former will depend on the magnitude of the ini- 
tial baryon number density ns(0) generated by CP- violating 
interactions at the bubble walls as well as particle number 
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changing reactions that extend into the unbroken phase where 
sphaleron transitions are unsuppressed. We will provide an 
expression for the BNPC that takes these considerations into 
account in Section VI. 

In organizing our discussion of these issues in the remain- 
der of the paper, we will be somewhat pedagogical since the 
issue of gauge-independence, though clearly of fundamental 
importance, is subtle. Consequently, some of the material 
consists of a review of earlier work, where we collect and re- 
cast in this context results that have been previously derived 
in the literature. In addition, we observe that although the 
complete evaluation of Fsph is vastly more difficult than de- 
termining Tc, it appears that at least in the SM the numerical 
impact of the treatment of Tc appears to be significantly more 
important than the one-loop contributions to Fsph associated 
with fluctuations about the sphaleron. Consequently, most of 
our discussion is focused on Tc and V^ff . 

To that end, in section II, we discuss the effective poten- 
tial and the Nielsen identities using a general model that ac- 
commodates an arbitrary extension of the SM. We identify 
the gauge-independent quantities derivable from the poten- 
tial. We also discuss how a failure to consistently imple- 
ment the Nielsen identities leads to gauge-dependence in Tc- 
We subsequently outline our gauge-independent treatment of 
the critical temperature in section III, and discuss the gauge- 
independent sphaleron rate in section IV. We provide a con- 
crete, numerical application of the gauge-independent analy- 
sis of the SM in section V, and we discuss the practical im- 
plications and theoretical limitations of our procedure - along 
with the modified BNPC - in section VI. In that section, we 
also examine the dependence of the BNPC on a variety of the- 
oretical inputs. Finally, we conclude in section VII. Several 
formal results are contained in the Appendix. 



II. GAUGE PROBLEM IN THE STANDARD ANALYSIS 

We begin our discussion by focusing on the effective poten- 
tial Veff, since it from this quantity that the critical temperature 
Tc and the background field (f>^i„ are derived. The standard 
practice has been to work in Landau gauge — 0), but we 
will leave ^ arbitrary to illustrate how gauge dependence of 
Tc and (f)^i„ appear in the standard analysis. For pedagogical 
purposes, we will first derive the zero-temperature effective 
potential in a general gauge and follow up with its general- 
ization to finite temperature. After briefly discussing implica- 
tions of the gauge problem in the standard analysis, we will 
discuss how Nielsen's identities shed light on this problem. 



sentation, the generators T", a = I, . . . , N are purely imag- 
inary and antisymmetric. For convenience a negative imag- 
inary unit, —i, is factored out, so that generators are purely 
real and antisymmetric T" = —i T^^. All gauge bosons of the 
theory are assembled into a vector, so that the generating 
functional and Lagrangian of the theory may be written in a 
compact manner: 



Z[j] = J p^PAe*-'''''^^^^'^') 



(8) 



Our convention for the gauge covariant derivative and field 
strength tensor is' L>^$i = 9^$^ + ^"(^r^J.y^j, and 
F^^ = d^A'i, - d^A^ ~ gf^^^A^^Al, where are the struc- 
ture constants of Q. The effective potential is obtained by 
passing from the generating functional of connected Greens 
functions W[i] = —ilnZ[j] to the generating functional of 
proper (IPI) vertices, F[(/)ci] via a Legendre transformation 



(10) 



where (pdix) = dW[j]/dj{x) is the classical field conju- 
gate to the source, j{x). For space-time homogenous clas- 
sical fields (pciix) = 0ch we obtain the effective potential, 
r('?^'ci) = ^(vol) Veff{(f>a)- To properly define the functional 
integral in (8), we must factor out a set of redundant gauge 
configurations. Doing so introduces the gauge-dependence 
that is the central focus of our analysis. 

To make this gauge-dependence explicit, we derive the ef- 
fective action following the background field method, where 
the scalar fields are (homogeneously) shifted as a change 
of variables under the functional integral in (8). We write 
^i{x) ~ (pi + (t>i{x), where (pi are space-time independent 
background fields^ and (pi [x) represent quantum fluctuations 
around the background field. 

Organized by powers of quantum fields, the resulting la- 
grangian after shifting and before gauge-fixing is 



C{x-3) = -V{^)+3^ + <P, 



1 



dV 



+ -0, [ - - 1/2(0)] 0, - a^0,A;i(gr>), 

- iA^[(9V'^ - df^dnS^" + m\{4>Y'g^^'^]Al, (11) 



A. Zero temperature effective potential 

We will work in a very general 4D theory[31] ("general 
model") containing an arbitrary number, n, of real scalar de- 
grees of freedom, which are assembled into a column vec- 
tor, ^i{x) (with i = 1, . . . ,n), transforming under a general 
(possibly reducible) real representation of an arbitrary (semi- 
simple) gauge group, Q with N generators. In the real repre- 



In general, different gauge coupling constants are assigned to each simple 
subgroup of Q. To avoid notational clutter, all coupling constants are iden- 
tically denoted as g and it should be clear which one we refer to, based on 
the generator or gauge field appearing adjacent to it. 
^ More precisely, one should write the background field as being dependent 
on the source, (f> = <p{j). The dependence on j is chosen so that, upon 
performing the Legendre transformation, 4> precisely coincides with (p^i, 
facilitating the evaluation of the perturbative loop expansion of the effective 
action r[</>,.i] [32]. 
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where 

MlU>)^^^V|^4>^^cj>,\;j, 

are field-dependent mass matrices for the scalar and vector 
gauge bosons, respectively. Terms cubic and higher order in 
quantum fields are not shown. 

We now impose the background-field dependent gauge con- 
dition^. 



(from now, we drop the 'cl' label and bars off the classical 
background fields): 



V;ff(0) = Kee + — r 

vol 



4- 



1 / dct0g,„(0) 

2 ^ Vdct 0231,(0) 



det0sc((/)) 
cletOsc(O) 
/ det0Fp(0) \ 



(17) 



where 



(12) 



on the quantum fields. This choice, motivated by the 
gauges, has the advantage that adding C^f = — (J-'°)'^/2^ to 
the Lagrangian cancels the —d^ (t>iA'^{gT°-(l))i mixing term 
appearing in (11). The compensating ghost Lagrangian is 



+ gr^id^vnv'A'^^ - a9T''<j>W'ri''{9T!,<Pj) • (13) 

Upon adding £gf and £gh to (11), the gauge-fixed generat- 
ing functional is 



Z\j]= j Vfj^VAV-qVii^e 



(14) 



with 



dV 



+ -</),[ -a2-il4(0) -^m^(0),,]</,,- 

+ ^^^[(9^5^'' - (1 - ^)d^d^)S''' + m\{4>T^g>^'']Al 

+ 'n'^''[-d^5'''' -^m\{^f'']if + , (15) 

where we have suggestively written 

^m\{^\, = ^{gT-^UgT'^^), (16) 

as the additional gauge-dependent scalar boson mass matrix 
arising from £gf . We prove in appendix B that since they are 
both built out of the same object, (gT°0i), the mass matrices 
m\{4))ij and m\{(j>)°''' share the same non-zero eigenvalues. 

To obtain the one-loop effective potential, we perform the 
Legendre transform, and the functional integrals are carried 
out in the gaussian approximation; terms cubic and higher or- 
der in quantum fields are dropped from (15). Formally, the 
result is a sum of functional logarithmic determinant ratios 



Ogau(0) = {d^g'^" - (1 - ^)d^d'')S-'' + m\{<i,Y<'g>''' 



(18) 



are the fluctuation differential operators appearing in (15). 
The functional determinants may be computed by going to 
Fourier space, and diagonalizing the mass matrices [35]. The 
result is 



Kff(0) = Kee(0) 



Trln(/_A/2.(0)_^„i^(0),^.) 



+ id-l) Trln (/ - mli^r") + Trln (/ - 



2 Trln {p^ - Cm\{(t)r^) - "free" 



(19) 



where the UV-divergent momentum integrals are dimension- 
ally regulated (d = 4 — 2e), and "free" refers to the field- 
independent subtractions arising from determinant ratios in 
(17). In our notation, the trace-log of a matrix means to 
sum over the logs of eigenvalues of the matrix. We note that 
each term in (19) represents fluctuations of specific degrees 
of freedom in the theory; the first term represents scalar and 
Goldstone fluctuations; the second arises from the transverse 
and longitudinal fluctuations of gauge fields, while the third 
term corresponds to the scalar (time-like) gauge field fluc- 
tuations; and the fourth term represents fluctuations of the 
Fadeev-Popov ghosts fields. 

After adding together the last two logs with identical argu- 
ments, we perform the momentum integrals using the standard 
formula (already expanded in powers of e). 



4(47r)' 



(2^ 



2\2 



ln(p2 



"fE — ln47r — In 



(20) 



^ It is more general to choose the fixed-vector gauge condition: dfj,A'^ " — 
^(piv!f = originally proposed in [33]. However, it was pointed out in 
[34] that the fixing condition does not allow for a consistent analysis of the 
vacuum structure under the homogeneity condition. 



The 1 /e poles are cancelled against appropriate counter-terms 
by choosing a renormalization scheme. We choose the conve- 
nient MS-scheme in which the 1/e — 7^ + In 47r are removed. 
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FIG. 1 : a. Gauge-dependent scalar-QED type graph describing (j)(j>- 
scattering that contributes to r[(^, A], b. External-leg corrections 
necessary to remove gauge-dependence from <j!)0-scattering ampli- 
tude, but absent from r[(;/), A]. 



The result is 



Kff(0) = H-ee(0) 
1 



E 

scalars.'i 



4(47r)2 



E 

gauge, a 



4(47r)^ 



In 



In 



m2 



E 

gauge, a 
n2 



4(47r) 



In 



(21) 



where m^i4>;(,) are eigenvalues of the scalar mass matrix 
Mfj (0) + ^ m\ ((/)) , and (0) are eigenvalues of the gauge- 
boson mass matrix m\{(j))°'^. The sums run over these eigen- 
values. 



and thus, the effective potential remains gauge-dependent. We 
show in Appendix B that in the specific case when (j) extrem- 
izes the tree level effective potential, the mass matrices decou- 
ple, leading to the necessary cancellation. 

The dependence of 14ff((/)) on ^ may be traced down to 
the dependence of the effective action, r[(j), A^J. We provide 
two equivalent ways of understanding the origin of its depen- 
dence on ^. First [32], the effective action, defined as the sum 
of 1-Pl graphs, includes a scalar QED-type graph shown in 
Fig. la, which is gauge dependent. The graph describes a 
physical scalar-scalar scattering process that cannot be gauge- 
dependent. The sort of graphs needed to cancel gauge depen- 
dence from the scattering amplitude is shown in Fig. lb. Since 
external leg corrections are not 1-Pl graphs, they do not enter 
the effective action, leading to gauge dependence. 

Second, the effective action is equivalently defined as the 
Legendre transformation of the generating functional of con- 
nected Green's functions W [j] with respect to the source j {x), 
where 



W[j] = -iln J V(f>VAVT]Vr] 



t j/d''x(£(e)+i,0,) 



(23) 



Unlike the gauge-fixed Lagrangian, C{^), the source term, 
(j>iji, appearing in the generating functional is not BRST in- 
variant since the source transforms as a gauge singlet. As a 
result, the effective action derived from the generating func- 
tional in this way will depend on the gauge-fixing procedure, 
and hence, the gauge parameter. 



B. Origin of gauge dependence 

Before generalizing to finite temperature, we pause to com- 
ment on the gauge-dependence of the zero-temperature effec- 
tive potential. The essential structure is shown in (19). There 
are three sources of the gauge parameter: from the scalar mass 
matrix S,rn\{(t))ij in the first log, from the time-like gauge 
boson mass matrix i_'mA{4>)°'^ in the third log, and from the 
Fadeev-Popov mass matrix ^771^1(0)"^ in the fourth log. As 
m\{(j))ij and (</>)"'' have identical non-vanishing eigen- 
values (appendix B), one could in principle arrange a cancel- 
lation among these logs by writing 

Trln (p2 _ A/2.(0) „ ^m^(0),,) 

Tr In (p2 _ A/2. (0)) + Xr In {p" - £,m\ (^),, ) (22) 

provided mj^{(p)ij and Mf-{4)) are simultaneously diagonal- 
izable and that their eigenvalues live in distinct subspaces. 
However, for general 0, the mass matrices do not decouple. 



C. Finite temperature effective potential 

We now apply these arguments to the finite temperature ef- 
fective potential. The minimizing field 4>mia{T) and the criti- 
cal temperature Tc are obtained from the finite temperature 
effective potential. The latter is derived from the partition 
function Z = Tr[e^^/^] that is expressed as a Euclidean path 
integral over field configurations satisfying periodic boundary 
conditions in Euclidean time. 

The generalization to finite temperature is achieved by re- 
peating the steps leading to (19), but with the energy integral 

replaced by a sum over Matusbara modes: / ^ T 

The potential is computed in the background field method 
with the appropriate gauge fixing as before. We mention 
that to retain BRST symmetry at finite temperature, the ghost 
fields must satisfy periodic boundary conditions, despite be- 
ing Grassman- valued [36]. After making these modifications, 
the generalization of (19) is 
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+ 00 

yeff(0,T)=Kee(<^) + - 



T 
2" 



(27r 



Trln {cjI + + Mf^{cf>) + ^m\[<P),,) 



+ {d-l) Trln (cj^ + p2 + m^(,^)«*') + Trln (cj^ + p2 + ^'nil'^)"'')) ^ 2 Trln (lj^ + p2 + ^m^(</))"'') - "free" 



(24) 



where w„ = {2TmT) are the bosonic Matsubara frequencies and where "free" is the free-field subtraction. After performing the 
frequency sum, the one-loop contribution separates into a zero-temperature part and a temperature-dependent part. 



scalar,?' 



, (25) 



where Vaei:{<l>) + Vcw(0) is the Coleman-Weinberg effective potential in Eq. (21), Jb{z'^) ~ /q°° dx ln(l — e v'^^+^^) is the 
bosonic thermal function, and where inf{(t); ^) and m'^{(f)) are, respectively, the eigenvalues of the scalar and gauge boson mass 
matrices as defined below Eq. (21). The thermal functions are typically numerically integrated or approximated using a Bessel 
function representation. 



Clearly, the finite-temperature effective potential (25) is 
gauge dependent. The dependence appears for precisely the 
same reasons described in the previous section. The scalar 
trace-log represented by the first term under the integral in 
Eq. (24) cannot be split in a manner similar to Eq. (22) be- 
cause the scalar mass matrices Mfj{(j)) and m\{(l))ij are gen- 
erally not simultaneously diagonalizable with eigenvalues liv- 
ing in distinct subspaces. 



D. Implications for the standard analysis 

We now show that the gauge-dependence of the effective 
potential leads to an unphysical treatment of the BNPC in the 
standard analysis. For the sake of clarity we temporarily re- 
strict our discussion to the SM. 

The critical temperature Tc marking a phase transition be- 



tween two phases 
by 



(1) 



and (f)'^' , derived from (25), is defined 



514ff 



Veff(c^7^c;0-v;ff( 



a(2) 



eif 



,Tc 



,Tc;O=0 (26) 
, = . (27) 

'■Tc 



In the SM, ip'^^^ = is the symmetric phase, while 0'^'' = 
4'min{Tc) characterizes the broken phase. In the standard anal- 
ysis, (26) and (27) are simultaneously solved to obtain both 
Tc and 4>mm{Tc)- Their ratio is compared to the bound (6) 
required for preservation of the baryon asymmetry. 

It is straightforward to see that naively inverting these equa- 
tions at one-loop order leads to a gauge-dependent estimate 
of the sphaleron rate. To that end, we consider the high- 
T approximation to the full effective potential. The thermal 
bosonic function Jb{z'^) admits a high-temperature expan- 
sion [26] 



-45 + 12^ 



2x3/2_ 



1 

32 



(28) 



with which the effective potential (25) may be cast as a poly- 
nomial in (h: 



(29) 



The coefficients D,Tq, E and A depend on the parameters of 
the underlying model. In the SM, the coefficients are [37] 



{29l + (gl + gi)'/') 



(30) 



967r 

and A = A + (^-dep. log) , 

where yt is the top yukawa coupling; gi and 172 are the U(l) 
and SU(2) l gauge coupling constants; and the scalar quartic 
self coupling A picks up a logarithmic ^-dependence. 

We observe that the coefficient of the quadratic term 
is gauge-independent, as one expects based on the gauge- 
independence of thermal masses (see e.g., Ref. [38]). In Ap- 
pendix C, we explicitly demonstrate this property for the gen- 
eral model. As we will discuss below, we take advantage of 
this property to define the high-temperature effective theory 
used to obtain a gauge-independent sphaleron scale. 

Unfortunately, the coefficient E is not only gauge- 
dependent but strongly so. For example, by choosing ^ = 3^/"^ 
the i?-coefficient can be made to vanish, and the barrier nec- 
essary for a first order phase transition is permanently ab- 
sent. One might hope that the ratio 4>mm{Tc)/Tc' removes 
this gauge-dependence. However it is straightforward to show 
that this hope is not realized. A simple calculation gives (for 
a pedagogical review, see Ref. [9]) 



Tc 



2E 
X 



(31) 



which remains gauge dependent. Hence the standard analysis 
leads to an unphysical treatment of the BNPC. 
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In various Standard Model extensions where either new 
scalar loops or tree-level operators generate large contribu- 
tions to the cubic term, it is conceivable that the impact of 
the gauge dependence is numerically minimized so long as 
one works in a region near ^ = 0. Nevertheless, such a sit- 
uation would be wholly unsatisfying as one would not know 
whether a perturbative computation near C = reflects reality 
or introduces a substantial gauge-dependent artifact. Instead, 
one should endeavor to maintain strict gauge-independence 
throughout the analysis in order to avoid any ambiguity. Be- 
fore moving on to our proposal for a gauge-independent anal- 
ysis, we briefly review some theoretical properties of the ef- 
fective potential that will prove useful later 

E. Nielsen identities 

The dependence of the effective action on the gauge-fixing 
parameter is described by the Nielsen identity [27] [34], and 
its precise form depends on the gauge fixing condition im- 
posed on the quantum fields. The identity follows directly 
from the BRST (non)invariance of the sourced generating 
functional. For the general class of linear gauges, the iden- 
tity reads 



Veil 




FIG. 2: A schematic illustration of the behavior of the exact effective 
potential as the gauge parameter ^ is varied according to Nielsen's 
identity. The values of the potential at its extrema stay unchanged 
but the fields extremizing the potential are gauge-dependent. 

gauge-independence of the critical temperature. However, as 
mentioned above, the precise value of the minimizing field 
(f>{Tc) remains gauge dependent (see Appendix E for details). 
Thus, the ratio (t){Tc)/Tc conventionally used to establish the 
BNPC is inherently gauge dependent. 

F. Artificial violations of Nielsen's identities 



^ = [ d'^xd'^y\Q{(l},A; x,y)-^^ 
J L ^4>^[x) 



(32) 



where, C and E stand for field-dependent vacuum correlators. 
The interpretation of the identity is clear: the value of the ef- 
fective action evaluated at its stationary points, i.e. at points 
where the classical fields satisfy 6r/S4ii{x) = (5r/(5yl° (.t) = 
0, is gauge-independent. This is as it should be since physical 
quantities are derived from the stationary points of the effec- 
tive action [39]. 

To compute the effective potential, the vector potentials are 
taken to vanish A^{x) = 0, and the scalar fields assumed to 
be homogenous (j>{x) = cj). In this case, the effective action 
reduces to the effective potential, r[0]/vol = — K:ff((/)), and 
the Nielsen identity simplifies to 



dV^ff dV,« 



(33) 



which carries the same interpretation: the effective potential is 
gauge-independent where it is stationary. Note, however, that 
the point in field space minimizing (or extremizing) the effec- 
tive potential 0min(C) is not gauge independent; as the gauge 
parameter is varied, the effective potential compresses or ex- 
pands along the 0-axis (see Fig. 2), while maintaining gauge 
independence of the energy density at its extremal points. 

The derivation of (32) and (33) does not rely on per- 
forming any non-trivial space-time or energy-momentum in- 
tegrations. Thus, Nielsen's identities are valid at finite tem- 
perature upon the standard replacement J ^ — !• Tj^n- 
An immediate consequence, proven in Appendix D, is the 



Although gauge independence of physical quantities is an 
exact statement following from Nielsen's identities, gauge de- 
pendence may be introduced to physical quantities as an arti- 
fact of an approximation scheme. In particular, naively trun- 
cating the perturbative effective potential at a finite order in 
perturbation theory leads to apparent violations of Nielsen's 
identity. In fact, we identify this naive truncation as the princi- 
pal source of gauge dependence in the standard determination 
ofTc. 

Such effects may be deduced directly from Nielsen's iden- 
tity. By expressing Vsff{(j)) and C{(j), £_) in (33) as a series in 

h, 

VM = Vo (0) + hVi (0) + V2{(b) + . . . (34) 

Ci(b,0 = Co + hcii(b) + re C2{<j)) + . . . , (35) 

and retaining terms through 0{h) in both sides of (33), we 
find 

dVo^.dVi dVo .( dVi dVo\ 

Since the tree-level potential Vb is strictly gauge independent, 
setting ©(fi,^) terms equal implies cq — 0. Upon setting 
0{h}^) terms equal we have 



dVi 



-Ci- 



dVo 



(37) 



that is, the one-loop potential is gauge-independent only 
where the tree-level potential is extremized, and not where the 
one-loop potential is extremized. Therefore, the critical tem- 
perature based on the one-loop extremum is gauge-dependent. 
We note that gauge dependence appears formally at a higher 
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order than the order of approximation. Yet, numerically, there 
is no limit to its sensitivity, and one should necessarily strive 
to obtain a critical temperature that strictly maintains gauge- 
independence at each order in perturbation theory. 



Here, (^i (T, ^) represents the genuine, albeit gauge- 
dependent, one-loop correction to the tree-level VEV (po at 
zero temperature, or to the background field at finite temper- 
ature. After substituting these back into (38) and expanding 
again in h, we find 



III. GAUGE-INDEPENDENT DETERMINATION OF Tc 

In the following sections, we detail one method for per- 
forming a gauge-independent perturbative analysis of the 
EWPT. As discussed above, exact expressions for the criti- 
cal temperature and sphaleron rate are not at our disposal; 
each has to be computed in an approximation scheme. We 
focus on the critical temperature in this section, and discuss 
the sphaleron scale in the next section. We stress that the 
computation of the critical temperature and the temperature- 
dependent sphaleron energy are two independent calculations, 
and that each should separately be gauge-independent. Since 
they are derived from different quantities - one from the effec- 
tive potential, the other from the effective action - the analysis 
should be treated as such, by computing both quantities sepa- 
rately. 



A. The h expansion 

The key to extracting a gauge-independent critical tempera- 
ture is to invert the defining equations (26) and (27) in a man- 
ner consistent with the ^.-expansion. To that end, we start by 
computing the effective potential to the desired order in per- 
turbation theory as outlined in Sec. II: 

Feff (0, T) = Voi(b) + hVi{(b,T) + h^V2i(b,T) + ... (38) 

We wish to minimize this function. 



We then write i 



as a series in h 



</'min = 00 + ft 01 (r, + ft' 02 (T, C) + 



(39) 



(40) 



where 0o represents any one of the minima of the tree-level 
effective potential. Upon substitution into (39), we find 



_dVi 



dVo 



ct>o+ri(t>i+. 

Ha0 



d'Vr 



0O + 'i01 + 

b 



902 



+ ...==0 
O{h^)^0, (41) 



where in the second line, we have expanded and organized 
in powers of H. The requirement that each coefficient must 
vanish individually determines the position of the minimum 
0niin at each order in perturbation theory; 



0(ft°) : 00 
0{h'): 0i(T,O 



m 

d(j) 



4>o 



d(j)'^ / <Po 90 



(42) 



^4ff(0min(T), T) = T/o(0o) + ftVl(0o,T) 

+ ft2 [1^2(00, T,O-i0?(T,O0Uo] +0{h^). (43) 



This formula for VI;ff(0niin(T), T) appears in Ref. [40], ap- 
plied to the high-temperature effective theory. 

At each order in h, the RHS of Eq. (43) is gauge indepen- 
dent in accordance with Nielsen's identity. An explicit ver- 
ification of this for the 0{h) term is as follows. The one- 
loop expression is shown in (24) and must be evaluated at the 
tree level minimum 0o. Thus, the mass matrices Mfj{(po) and 
TO^(0o)y in (24) are evaluated at 0o. At that point, these 
mass matrices are simultaneously diagonalizable and eigen- 
values live in distinct subspaces (Appendix B). As discussed 
in section II B, this situation allows us to split the scalar log 
into two terms as in (22) and finally achieve the needed can- 
cellation among the gauge-dependent logarithms. 

Equation (43) provides a gauge-independent closed-form 
expression for the value of the effective potential at its local 
minimum (or maximum) as a function of temperature. Since 
the function may be computed from any extremum 0o derived 
from the tree-level potential Vq, equation (43) defines a family 
of curves which follow the temperature evolution of the ex- 
tremal points. At each intersection point, the degeneracy con- 
dition (26) is satisfied, and that point potentially corresponds 
to a critical temperature. Since phase transitions typically take 
the system from one global minimum of Vgff to another, criti- 
cal temperatures occur when the two lowest curves intersect. 
Fig. 3 depicts a scenario for a hypothetical theory with three 
phases. Evolution of the universe proceeds from right to left. 
In this example, the phase transition occurs in two steps. First, 
a transition occurs from phase 1 to phase 2 at temperature 
Tc 1 . Then a second transition takes place from phase 2 to 
phase 3 at temperature Tc,2- 

Note that, in stark contrast with the current analysis which 
focuses on the complicated multivariable function 14ff(0, T) 
to follow the free energy of each phase, the foregoing method 
deals with a family of functions dependent on the single vari- 
able T. Because numerical algorithms that find intersection 
points of two curves are significantly more efficient than those 
that find global minima of multivariable functions, our method 
to find critical temperatures is substantially faster than the 
standard methods. This is an added benefit of the gauge- 
independent approach discussed above, as it would speed up 
the analysis of multi-step phase transitions in complicated ex- 
tensions of the SM that are subject to extensive parameter 
scans. 
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Vcff(«4min) 




T 



FIG. 3: Value of Kff(0rain) for three different phases as a function 
of T as determined by (43) in a hypothetical theory. Evolution of 
the universe proceeds from right to left, following the direction of 
the arrows. Intersections of the two lowest curves define the critical 
temperatures (there are two in this example). 

B. Qualitative analysis of the gauge-independent critical 
temperature 

One may get a sense of the behavior of the critical tem- 
perature at 0{h) by analyzing (43) in the high-T limit. We 
shall restrict ourselves to a single degree of freedom and study 
a phase transition from the symmetric phase to the broken 
phase. We normalize the potential such that it vanishes at the 
origin at all temperatures, so that the degeneracy condition 
(23) to 0{h) reads 

V,s{4>, Tc) = Vq{4>o) + ftl/i(0o, T) = . (44) 

After applying the high-temperature expansion of the thermal 
bosonic function (28) to the LHS of (44) we have 

D{T^ - T^)<i>l - ET4>1 - - . (45) 

Note that gauge-independence of this expression is ensured 
only when the fields are evaluated at their tree-level minima 
00 = Tf)^/2D/\ as dictated by the ^-expansion. Solving 
this equation for T yields two roots; only one gives a positive 
critical temperature 




where D was eliminated in favor of 0o- Assuming no tree- 
level cubic terms, we may take E small, since it is loop- 
suppressed. We the find 

rc = ^«0.71To, (47) 

that the critical temperature is approximately 30% lower than 
To (the critical temperature in the second-order limit). Our 
one-loop numerical study of Tq presented in Section V con- 
firms these expectations. Writing To = (f)a\/\/2D, we find 
the scaling behavior of the critical temperature 

Tc^fvf- ^^^^ 



As the quartic coupling A is increased, the critical temperature 
rises, too. 



C. Higher-order contributions 

When analyzing the behavior of the effective potential, it is 
instructive to investigate the impact of higher-order contribu- 
tions and to explore the limits of validity of perturbation the- 
ory. In this section we study these higher-order contributions 
to assess their importance to (43). 

It is well-known that at high temperatures the perturbative 
expansion breaks down due to uncontrolled g^T^/m^ terms 
appearing at higher orders[26][41][42]. These "plasma damp- 
ing" contributions are induced by fluctuations of ?? > Mat- 
subara modes at high temperature. This malady is cured by re- 
summing these terms to all orders to obtain a "ring-improved" 
effective potential V^™^ (</>, T) that includes the damping ef- 
fects. In this subsection, we will discuss two possible gauge- 
independent methods of including such effects. We will then 
conclude by briefly discussing the 0{h^) term of (43). 



1. High-T effective theory 

One may include damping effects by going to the high-T ef- 
fective theory by integrating out the heavy Matsubara modes 
via dimensional reduction. Nielsen's identities (32) and (33) 
apply to any gauge theory, including effective theories. Pro- 
vided no explicit gauge-dependence is introduced in the Wil- 
son coefficients of the effective theory, stationary points of 
the effective potential derived from it will remain gauge- 
independent. For the same reasons that the effective action 
is gauge-dependent, so is the effective potential obtained with 
the dimensional reduction procedure. Indeed, the Wilson co- 
efficients of all terms, apart from the C(T^) terms are gauge- 
dependent. However, by working in the limit T rriy^^, one 
can justify retaining only these 0{T'^) terms of the effective 
theory. These terms are gauge independent because thermal 
masses are gauge independent. In principle, then, it is possible 
within this limit to derive a gauge-independent critical tem- 
perature from the high-T effective theory with the g^T^ jrr? 
terms implicitly re-summed. 

This procedure has been applied to the SU(2) Higgs model 
in Ref. [40] , where it was argued that the fi expansion for Tq 
breaks down at two-loop order. We revisit this analysis briefly 
below and also conclude that this approach is not useful for 
obtaining information about the phase transition. 

To proceed, we start with the effective lagrangian, where 
the heavy Matsubara modes are integrated out[40, 43, 44]: 

^3 = ^0<i>.f + m + ^5],,(T) 

+ k'"^"'^ AIA\A%A% + iJf/AgA(5$,$j- + , (49) 
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where F is the 3D field strength tensor, I](r) and n(T) are 
scalar and gauge-boson thermal mass matrices, respectively; 

gives the tree-level quartic coupling of the scalar and 
time-like gauge fields; and A"'"^'' is a T-independent loop in- 
duced quartic interaction of the time-like gauge fields. Note 
that we have not yet performed the integration over the 71 = 
Matsubara modes. In this theory, the tree-level scalar potential 



V 



high-T 







($,T) = V{^) 



(50) 



In the standard SU(2) model with a single degree of freedom, 
this potential may be written in a form similar to the high 
temperature-expansion of the effective potential derived from 
the full theory (29): 



V 



high-T ( 







f), T) = D{T'^ - T^)(j)^ 



(51) 



where DTq ~ — /^^/2 and where DT^ incorporates the effect 
of the scalar thermal mass with D being given in (30). In 
contrast to the high-T expansion of the full gauge-dependent 
effective potential (29), the RHS of Eq. (51) contains no T^^ 
term. Its appearance in (29) is due to Matsubara zero modes, 
which we have not yet integrated out to derive (51). Thus, at 
tree level, the phase transition is second order and the critical 
temperature is Tq. 

The one-loop effective potential is computed following the 
background field method, with an Rc gauge fixing condition 
similar to (12). Upon performing the functional integration 
over the zero modes, the high-T effective potential acquires 
the non-analytic terms necessary for a first order phase transi- 
tion. 



high-T/ 



eff V/^) T) — 



_ Trhigh-T 



T 

12^ 



E 



3/2 



e-values 

2 /j,\a6\3/2 I /2 /r\ab , Tiab /rr\\3/2 



ah\3/2 



(52) 



where the sum runs over the eigenvalues of the field- 
dependent mass matrices defined in below (11) and (15). 

The potential is gauge dependent as expected. Neverthe- 
less, one might expect to obtain a gauge independent critical 
temperature by applying the ?i-expansion (43) on the high- 
T effective potential and determining intersection points as 
outlined section III A. Indeed, the cancellation of gauge- 
dependence operates in a manner analogous to that in the 
full 4D theory. The difference in this case is the tree-level 
extrema derived from (50) are now temperature dependent 
(j>o = 4'q{T). So, when the scalar mass matrices are evalu- 
ated at these points in accordance with (43), the eigenvalues 
of Mf^{MT)) + ^^AT) and ^m\{<j>^{T)) decouple, al- 



lowing us to write 

[M^.(0o(T)) +^mi(0o(T))^, + E.,(T)] 



3/2 



i3/2 



(53) 



and arrange a cancellation with the last term of (52) after ap- 
plying theorem 2 in Appendix B. 

Before discussing the physical implications of the above 
procedure, we contrast this method with the approach taken 
in Ref. [40], which, in addition to V^ff and 4)mia, Tq is itself 
written in powers of h: 



Tc=n + hTi + h^T2 + 



(54) 



As before, the leading order critical temperature is To and 
the transition is second order Inclusion of the 0{h) terms 
yields a contribution Ti proportional to the product"* of H"^' 
and n"''(T). At 0{It?), the scalar mass acquires a logarith- 
mic dependence on the renormalization scale yu as needed to 
maintain the renormalization group invariance of the poten- 
tial. Ref. [40] finds that this is associated with a two-loop 
contribution to the potential of the form 



In 



m2 



(55) 



When solving for T2, the argument of the logarithm must be 
evaluated at T = To in order to maintain consistency with 
the ft-expansion of Tc- In this case, the argument of the log- 
arithm vanishes, signaling the presence of an IR divergence. 
As Ref. [40] concludes, the h expansion clearly fails when 
implemented in this manner. 

We believe that our method proposed in section III A does 
not suffer from this breakdown. Note that our method of de- 
termining the critical temperature just requires the analysis of 
the h series expansion of Kff ((/)min) in (43). Since this quantity 
is already gauge-independent at any value of the temperature, 
a series expansion of Tc is not needed. Instead we take the 
temperature as an external parameter and then numerically 
solve for Tc that satisfies the degeneracy condition (26). In 
general, then, the value of Tc appearing in the argument of 
the logarithm in Eq. (55) will differ from its leading order 
value To, and so one should not encounter an IR divergence 
with this method. One must of course separately treat the spe- 
cific case when the solution of the degeneracy condition gives 
Tc = To exactly. As we also note below, we do not encounter 
any IR divergences at 0{fP') when determining Tc using the 
full theory in the high-T regime. Thus, it appears to us that 
the breakdown of the h expansion as obtained in Ref. [40] is a 
consequence of expanding Tc itself. 

We emphasize that these observations do not constitute a 
criticism of the work of Ref. [40] . In that work, the full 0{h?) 
contributions to the minima of the effective potential in the di- 
mensionally reduced effective theory were computed explic- 
itly and shown to be gauge-independent and an expression for 



' In the notation of [40] these are denoted and m^) , respectively. 
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Q' D 



FIG. 4: A representative ring diagram contributing to AVnng{4', T)- 
Tlie Matsubara n = Q mode propagates in tiie lower loop. The 
smaller loops are hard thermal loops, each of which contributes a 
factor of g^r^/m^. 



the value that minimum similar to Eq. (43) appears. Moreover, 
a careful discussion of IR divergences for T in the vicinity of 
To was given. We merely differ on the application of the h- 
expansion to Tc itself. 

That being said, we also conclude that this approach using 
the dimensionally-reduced, high-T effective theory is unlikely 
to yield useful information about the true critical temperature. 
The reason is that in the dimensionally reduced theory, at lead- 
ing order in h in the SM, there exists only one minimum at 
any given temperature. At T > Tq, this minimum lies at 
(^0=0. For T < To, (/>o = gives a local maximum while 
= 2D{T^ - T^)/X. At T = To the two coincide, and 
one encounters the IR divergence when evaluating the higher 
order terms at (po ~ 0. It is unlikely that for T < the 
higher order corrections convert the tree-level local maximum 
at the origin into a local minimum, so that one expects the 
only critical temperature that ever appears is To. Thus, we 
abandon this approach for incorporating the plasma damping 
effects into the determination of Tc. 



2. A ring-sum prescription 

The second alternative to including the damping effects en- 
tails a variant on the conventional approach based on the full 
4D theory. In this approach, the g^T^/m^ terms are under- 
stood to arise from the high temperature limit of ring (or daisy) 
graphs shown in Fig. 4. The summation of these graphs to all 
orders results in an additional contribution to the potential (25) 
of the form 



,T) = V^M<P,T)~V^^{^,T), 



(56) 



where 



e-values 



3/2 



(57) 



is the non-analytic part of Kff T) in (C4) and gets replaced 
by the thermal damping expression 



(A4(0)+C™A(0k+S.,(T)) 

e-values 

+ 2(m^((/.)'^'')'/' + (mi(0)"^ + n"^(T))'/' 



3/2 



ab\3/2 



(58) 



which coincides with the 1-loop expression of V^g^^'^{4>, T) 



in (52). These may be added to bring Aiding ( 
recognizable form 



6, T) into a more 



nngV 



6,T) = 



- h 



T 
'12^ 



scalars, i 



(59) 



longit. 
gauge, a 



where m'^{(j); ^) — A'P{<j>) + £^mA{4>) is the scalar mass ma- 
trix, mA{(t>) is the gauge boson mass matrix, and S(T) and 
n(T) are their gauge-independent thermal self-energy matri- 
ces which start at 0{h). The sums run over the eigenvalues of 
these matrices. 

Note that since S(r) and n(r) are 0{h), the procedure of 
re-summing these graphs necessarily departs from the strict h- 
expansion needed to ensure gauge-independence. Thus, it is 
a non-trivial problem to implement the re-summation into the 
gauge-independent analysis as one would have to analyze the 
all-orders expression of (43) and identify the g^T^ /m? terms 
that need to be re-summed. 

Instead we provide a gauge-independent prescription for 
its inclusion based on the following observation. Earlier, we 
pointed out that vj^^^ (0, T) is just the one-loop expression 
of the effective potential derived from the high-T theory, and 
V\^^ (0, T) is the 0(T) part of the effective potential derived 
from the full 4D theory. Thus, if we evaluate V^^J at (po{T), 

which minimizes the tree-level high-T potential V^^ (<^, T) 
in (50), we may apply (53) to cancel the gauge-dependence in 



[B], 



)(T),T). Similarly, if we evaluate vj^^ at which 
minimizes the tree-level potential of the full 4D theory (9), we 
may apply 

(A4(</>o)+^m^(</>o).,)'^'-^ 

the 4D analog of (53), to cancel the gauge-dependence in 
i>o,T). Therefore, the prescription we propose for ob- 



lA], 



taining the gauge-independent expression for AVnng is 



AV:-^iT)^vj^iMnT)~K, 



[A], 
ring ^ 
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so that the expression for the gauge-independent one-loop 
ring-improved effective potential at its minima is 

V-eff (0™n (T) , T) = (<^o ) + ^^1 ((^0 , T) + AV,f,^ (T) .(61) 

As it stands, the construction of AV^^^l' is somewhat ad hoc 
as we have not yet rigorously demonstrated that (60) is the 
result of re-summing the g^T^/m'^ terms in (43) to all orders 
in h. Instead, using a simple 0'* theory, we illustrate that our 
prescription correctly gives the hard thermal loop contribution 
to the 0{h^) term of the series in (43). 

Starting with the tree-level potential of 0'* theory. 



A 



2^ 4^ 
we find that the tree-level minimum 
dependent mass m'^{(j)) are given by 

d'Vo 



00 = vV/A 



(62) 

00 and the field- 
/i^ + 3A02 (63) 



Following the procedure outlined in section II, a little work 
shows that the one-loop finite temperature effective potential 
(25) is given by 

Vi{<P,T) = l/cw(0) + ^JB{m\(^)/T^) 



90 



J.2 



T 



24 127r 



(64) 



from which we derive the one-loop correction to the VEV 
(42), 



8fi 



3V2A 
Sir 



T + 0{T°) 



(65) 



We also find the tree-level potential (50) of the high tem- 
perature effective theory to be 

,r) = i(-/i2 + S(T))02 + ^04^ (66) 

where — hXT^/A is the thermal mass. The tree-level 

minimum of this potential is given by 



^high-T^ 



0o(T) = y'(/.2-i](r))/A. 



(67) 



Then, the result of re-summing the ring diagrams Fig. 4 gives 
(56) 



nngW, T) — V^j 

T 



ring ^ 



-h- 



[A], 
ring ^ 

3/2 



,T) 



3/2 



(68) 



The prescription to construct the gauge-independent poten- 
tial is to evaluate the first term term in (68) at (pQ (T) given by 
(67) and the second term at (po given by (63). The result is 



nng 



^12^ 



2^2 



XT' 



3/2 



^2 V y3 

8V27r 



(2^2)3/2 

(69) 



We now show that this 0{H'^T^) term is precisely what is con- 
tained in the 0{h^) term of (43) 



(70) 



An explicit calculation of the two-loop bubble graph of the 
type in Fig. 4 contributing to V2(0o) yields 



V2(0O,2")|bubb. 



O 



-V_y3 



(71) 



Then, substituting (63), (65) and (71) into (70) we find 
(eqn.70)^... + ^T3 + (.r), 



(72) 



that the 0{h^T^) term is in agreement with (69). 

We have carried out a similar analysis, drawing on the re- 
sults in Ref. [45], to show that the prescription given in (60) 
yields the 0{h?T^) in the SM. We remind the reader that 
while this result is not a rigorous demonstration of the va- 
lidity of this prescription, it does indicate that the gauge- 
independent quantity AVj^|^jJ gives an all-orders estimate of 
the effects of plasma screening which correctly reproduces the 
corresponding 0{h^T^) contribution arising from the full po- 
tential. In section V A, we will study the numerical impact of 
including AV^^^- when determining Tc- 



3. Two-loop order 

In addition to investigating the all-orders re-summation of 
the plasma damping corrections, one may also inquire about 
the remaining higher-order corrections. To that end, we con- 
sider the explicit 0{h^) terms in Eq. (43). We first observe 
that this contribution contains the term 



5^1/0, 



(73) 



which is manifestly negative due to the concavity of Vb at the 
minimum. By itself, this negative contribution will lead to 
an increase in Tc, and is therefore suggestive of the impor- 
tance of remaining 0{ti?) contributions. Of course, gauge- 
invariance requires inclusion of the explicit two loop term 
V2{(t)o,T). At present, we are not aware of any computation 
of V2 in an arbitrary gauge in the full theory^ and are unable 
to rigorously demonstrate the gauge-invariance determination 
of the 0{K^) contribution. 

On the other hand, the authors of Ref. [45] have computed 
this contribution in the full theory for the SM using the Lan- 
dau gauge. Thus, use of their results in combination with a 



^ The corresponding computation in the dimensionally-reduced, high-T ef- 
fective theory has been canied out in Ref. [40]; see above. 
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computation of (73) in the SM with same gauge should yield a 
gauge-independent result. Unfortunately, only expressions for 
V2 {(j), T) in the high-T expansion are available from Ref. [45], 
so it is not possible in the present study to provide a gauge- 
independent determination of the complete 0{h^) contribu- 
tions. We again defer a full two-loop computation to future 
study. Nevertheless, to the extent that the high-T approxi- 
mation to V2{(t),T) provides a numerically realistic estimate 
of the full result, one may use the former - together with our 
computation of (73) to obtain an indication of complete 0{h^) 
impact. In Section V A we show that the trend from the 0{fP) 
terms is to raise Tc, though the precise value should be taken 
with a grain of salt. 



IV. GAUGE-INDEPENDENT SPHALERON ENERGY 
SCALE 



The second factor appearing in the BNPC (6) is the 
sphaleron scale, commonly taken to be the scalar field 
4>mm{Tc') minimizing the effective potential. As empha- 
sized above, Nielsen's identity for 14ff (0, T) implies that the 
background field is an inherently gauge-dependent quantity. 
Therefore, we abandon (j)min{Tc), and we revisit the compu- 
tation of the sphaleron rate in search for a gauge-independent 
scale. In doing so, we recapitulate the work of Refs. [23, 
46, 47], wherein an appropriately gauge-independent formu- 
lation of a perturbative sphaleron rate computation is devel- 
oped. Along the way, we point out what we believe has mo- 
tivated a departure from this framework, leading to the intro- 
duction of additional gauge-dependence via the appearance 
of <^min(7c) in the BNPC. We then advocate a return to the 
gauge-independent treatment. 

The sphaleron rate Fsph is proportional to the free energy 
of the sphaleron gas g . For a dilute gas, it is related to 
the effective action for one sphaleron Feff [(^'^p^] normalized by 
the energy density of the electroweak broken phase Feff [0^^] 
[46]: 

Tsph = — Im Fs.c,. = — Im — !^ 



Im e-rcrrl^^P'Vrcrrl./.™] 



(74) 



where a;_ is the frequency associated with the unstable mode 
of the sphaleron. Since the sphaleron represents a saddle point 
of the energy functional, the rate is expected to be gauge- 
independent on account of Nielsen's identity. In order to en- 
sure gauge-independence within the perturbative context, the 
effective action must be computed in manner that is consistent 
with the ft-expansion in analogy with (43): 



A. Dimensional reduction 

For any extension of the SM, the electroweak sphaleron 
lives in an SU(2) subgroup of the full gauge group. Hence, we 
shall simplify our analysis by restricting ourselves to the mini- 
mal standard SU(2) model. The task at hand is to compute the 
path integral in (8) and (10) around the sphaleron background. 
In general, including finite temperature loop contributions to 
Tefti^^''''] is considerably more challenging than doing so for 
Vcs{<f>, T). The calculation may be simplified, however, when 
temperatures are sufficiently high such that the n > {) Mat- 
subara modes may safely be integrated out. As discussed in 
the previous section, this dimensional reduction procedure is 
gauge-dependent, but the Wilson coefficients of the O(T^) 
terms are gauge-independent. Gauge independence is main- 
tained by retaining just these terms in the effective theory. 

After imposing the temporal-axial gauge W^=o = on the 
SU(2) gauge fields W^, the path integral for static field con- 
figurations reads 



^3 

with 



\D,H\ 



1, 



W.t,Wtj+ViH,T) 



(76) 
(77) 



V{H, T) = 2D{T'^ - T^)H^H + \{H^ Hf 
^ X{H^H -v{Tf/2f , 



(78) 



where as before 2DT^ is the Higgs thermal mass; 2DTq = 
jjp' is the zero-temperature Higgs mass parameter; and 



v{T) = v,Jl-T^/T^ 



(79) 



with = 246 GeV being the tree-level expectation value of 
the Higgs field. Notice that H'^ = (0, v{T)/V2) mmimizes 
V{H, T) for T < To. It is over this action that the path inte- 
gral is performed to obtain the sphaleron rate. 



B. Tree-level sphaleron energy AiJsph 

We search for the tree-level sphaleron solution around 
which the path integral is to be computed. The Klinkham- 
mer and Manton[48] ansatz for the sphaleron, motivated by 
the electroweak instanton, is^ 



g X' 



(75) 



v{T)/V2 



(80) 
(81) 



Here, Slcj)"^ ] is the tree-level action and ri[0g'' ] is the one- 
loop fluctuation determinant, both evaluated around the tree- 
level sphaleron configuration ^q''''. 



* The minus sign in the ansatz for is due to our sign convention of the 
covariant derivative Di = di + igWt. 
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which we collectively denote as (f/^f^. Here, /(x^) and /i(x^) 
are the gauge and Higgs field radial profile functions, respec- 
tively. They are determined by imposing the stationary condi- 
tions on the energy functional and satisfy the boundary condi- 
tions, 



/(o) = Mo) = o 

/(oo) = Moo) = l 



(82) 
(83) 



Substituting (80) and (81) into (77) yields the energy func- 
tional in terms of the radial profile functions. 



ph 



T 



(84) 



with 



f2 /d/l\2 



(85) 



where r = gv{T)\:x.\ is a dimensionless radial coordinate. The 
stationary condition for the radial profile functions are 



dr 



df^ 
, dh 



dr 



2/(1 2/) 
A 



1 



f'h'il^f) (86) 



9^ 



(87) 



and must be solved numerically. With the radial profile 
functions, /(f) and h{f) known, the function B{X/g^) may 
be numerically evaluated, providing a quantitative result for 
S'3[(/)g''^]. The dimensionless function B{\/g^) is an 0(1) 
weakly varying function of its argument [48]. Therefore, the 
quantity which sets the scale of the sphaleron energy is v{T). 

Two points about v{T) and Tq should be emphasized. First, 
because the scale v{T) introduced by dimensional reduc- 
tion minimizes V{H,T), it is tempting to interpret v{T) as 
4'min{T), the value background Higgs field at temperature T. 
As discussed above, however, they only coincide in the high- 
T effective theory where the T(p'^ terms essential to bubble 
nucleation are absent. Consequently, it is best within the con- 
text of the calculation of the sphaleron rate to treat v{T) as an 
auxiliary quantity that sets the scale of the sphaleron energy. 

Similarly, the temperature Tq that appears in the expres- 
sion (79) for v{T) is not the same as the critical tempera- 
ture. Rather, Tq is the temperature at which the sphaleron en- 
ergy (84) vanishes at tree-level in the dimensionally reduced, 
high-T effective theory. At sufficiently high temperatures 
(T > Tq), the sphaleron rate contains no exponential sup- 
pression that is crucial for preservation of the baryon asym- 
metry. As the transition becomes more weakly first order Tc 
approaches Tq from below, leading to greater washout of the 
baryon asymmetry. 



C. One-loop fluctuation determinant 



the background field is the tree-level sphaleron configuration 
(f)'^^ determined above. The functional integration is per- 
formed in the gaussian approximation, formally giving ratios 
of functional determinants 



- eff 



J,sphl 



S3 



^1 



detao,((/)f ) 
detObos(</''^*) 



hill 



dct0Fp(0f ) 

detOFp(,^EW) 



(88) 



where Obos and Opp are the bosonic and ghost fluctuation op- 
erators, whose eigenvalues are proportional to {v/gTy^'^ as a 
result of coordinate rescaling. 

A precise calculation of these determinants is challenging, 
but the essential structure is as follows. The sphaleron back- 
ground breaks three rotational and three translational symme- 
tries of the system, giving rise to six zero modes in the spec- 
trum of Obos(0o'^)- The integration over these modes give rise 
to volume pref actors (AfV) and six factors of the proportional- 
ity constant {v /gT)^/'^. Here, M and V are normalization and 
dimensionless integration volume factors, respectively. Let- 
ting K denote the ratio of fluctuation determinants with zero 
modes removed, the sphaleron rate (74) then takes the form 



sph 



27r 



-A£;,pi,/T 



UJ , ,3/' \ „.-AE,ph/T 



= i^WV\{NVU{gvyy— I Ke 



(89) 



where in the second line we have expressed the dimensionless 
volume factor Vtr in terms of the physical volume using 
f = gv{T)\x.\. The factor of 1/2 arises from an analytic con- 
tinuation of the integration over the negative mode. 

The technically involved and numerically intensive task of 
evaluating k for the standard SU(2) model was achieved by 
Carson, Li, McLerran and Wang[23], and by Baacke and 
Junker[24, 25]. To our knowledge, similar perturbative com- 
putations do not exist for SM extensions. Rather, the conven- 
tional practice has been to attempt to approximate the one- 
loop contribution by replacing the potential V{H, T) in (78) 
with the one-loop effective potential Kff (0, T) in (25) derived 
from the full theory. In principle, doing so includes some 
part of the one-loop fluctuation determinants, but clearly ne- 
glects contributions involving spatial variations in the scalar 
and gauge fields. 

Following this practice and repeating the steps leading to 
(84) gives a similar result but with the gauge-independent 
scale v{T) replaced by the gauge-dependent Higgs field 
4>mmiT), leading to 



r.ph ^ ^(AAy),,(AAV),o,(5<^mi„)= , 
where 



The calculation of one-loop effective action F 



eff I 



isph 



)is 

performed following the background field method, in which 



T ^ gT 
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For a given choice of gauge, 4>min{T) naively approximates 
the full background field more closely than does v{T), so it 
may seem plausible that the replacement would lead to a bet- 
ter estimation of the sphaleron rate. By integrating the rate 
law using this expression for the sphaleron rate to derive the 
washout criteria (details in Appendix A), one arrives at the 
criterion in equation (6). 

In view of our discussion of T4ff (0, T) it is immediately 
apparent that making the foregoing replacement introduces 
gauge-dependence into AE^ph/T — and thereby into Fsph — 
at any temperature. Moreover, as observed in Ref. [49] in 
the context of false vacuum decay, this procedure double 
counts the n = Matsubara contributions to the sphaleron 
rate. Since these degrees of freedom already contribute 
non-analytic terms to Vcff {(t),T) from which (pmmiT) is ob- 
tained, including them in the fluctuation determinant calcula- 
tion would count them twice. 

In light of these two problems, we advocate abandoning the 
ad hoc replacement V{H, T) — > Kff (0, T) in favor of adher- 
ing to the gauge-independent formulation of Ref. [46] and 
using v{T) instead of (i)m\n{T) to characterize the sphaleron 
rate. When attempting to approximately establish the effi- 
cacy of baryon number preservation, one should then use the 
gauge-independent ratio of v{Tc)/Tc with Tc computed us- 
ing the /i-expansion as above and v{T) obtained from the 
dimensionally-reduced effective theory in the high-T approx- 
imation. Including 0{h) contributions to Fsph associated with 
Matsubara zero mode fluctuations about the sphaleron the re- 
quires a careful evaluation of the functional determinant. 



V. AN ILLUSTRATIVE APPLICATION: v{Tc)/Tc IN THE 
STANDARD MODEL 

As an illustrative application, we now apply our gauge- 
independent method to analyze the EWPT at 0{h) in the SM, 
and estimate the impact of higher order contributions. Since 
there exist Monte Carlo results for both Tc and Fsph in the 
SM, we may compare the corresponding quantities obtained 
with the gauge-independent perturbative analysis to assess the 
reliability of the latter As a result, we will conclude that in- 
clusion of two-loop contributions to Kff (0, T) is likely to be 
essential for obtaining a reasonable value of Tc- Later, when 
discussing the criteria for preservation of the baryon asym- 
metry (section VI), we will argue that a reliable perturbative 
computation of the fluctuation determinant is also likely to be 
important, though probably less decisive than the evaluation 
of Tc- We also argue that even when gauge-independence is 
maintained, a perturbative treatment is at best indicative and 
that definitive statements about baryon number preservation 
are likely to require non-perturbative studies. 

The two ingredients needed for this analysis are the crit- 
ical temperature and the sphaleron scale. While contribu- 
tions from fermions have been previously omitted, our nu- 
merical results include effects from the top quark. We be- 
gin by computing the effective potential for the neutral Higgs 
field, including the effects of the top quark. In the standard 



SU(2) X U(l) model, the four real degrees of freedom 



H 



1 f ^1+ 



(90) 



are placed into a vector $j = $2, ^3, ^i) - Explicitly, 
the generators T^g in the real representation {—i factored out) 
are 



rp3 _ 1 

-'Re ^ 2 



Tl = 



^Re — 2 



-1 
-1 



where the {T^^, T^^, T^^} are generators of isospin and T^^ 
is hypercharge. The Euclidean Lagrangian takes the form 



1 9 



with 



(91) 
(92) 



where D^^d^+ gT^^W^, and - {W^,,W'^, Wj^, B^}- 
We shift the scalar fields following the background field 
method. 



^ = (0, 0, h, 0) , 



(93) 
(94) 



and we immediately obtain the tree-level potential in terms of 
the classical field h, 



Voih) 



2" 4 

We record here the extrema of the tree-level potential: 



(95) 



h, 



(2) _ 



.0 - ± VmV^ = ±246 GeV. (96) 
Next we impose the gauge-fixing condition, d^^W^ + 
£.(f>i{gT^s4')i = on the quantum fields. Retaining terms 
quadratic in fluctuations, we read off the field-dependent mass 
matrices necessary to construct the effective potential; 

■2 ' Xh^ 



3Xh^ 



-^2 + Xh^ 



{9" 



(9' 



\ 



9 

-99' 



-99 

„'2 



(97) 



(98) 



(99) 



After abbreviating their eigenvalues 



TOq(/i) = -^^ + A/i^ , m^{h)^\g'^h^ 



ml{h) 



0. 



(100) 
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we write down the 1-loop finite-temperature effective poten- 
tial (25) 



A=0.035 {mM~65 GeV) 



(101) 



where explicit expressions in the MS renormalization scheme 
are provided in Appendix F. To determine the critical tem- 
perature, we use (43) to follow the temperature-evolution of 



the free energy of the symmetric /ig^-* and broken h]f' phases. 
Then the 0{h) degeneracy condition defining the critical tem- 
perature of the phase transition reads 



(2) 



Tc). 
(102) 



Before we proceed with numerical results, we explicitly 
show that both sides of (102) are gauge-independent, refer- 
ring to (Fl) and (F2). The statement is trivial for the LHS, 
since the one-loop potentials V^^^{h) and V^^'^ {h,T) are 

evaluated at the origin /iq^"* = 0. At that point, and 7ti^ 
both vanish and all gauge-dependent terms that multiply them 
disappear. In the RHS, the potential is evaluated at the tree- 
level broken phase minimum /iq^' — y^JjFJX. In this case 
rriQ = 0, and the gauge-dependent terms cancel. 

Let us illustrate how theorem 1 in Appendix B operates 
for the matrices in eqns. (97) and (98). Notice that the 
sum Mfj{h) + S^m\{h)ij that goes into the computation of 
the effective potential (24), for general values of the Higgs 
background field h, has mixed dependence on the gauge- 
parameter. This blocks us from being able to split logarithms 
as in (22). However when evaluated at the tree-level minimum 
= ±v7iV^, we find 

Mf^.(4'') = diag(0, 0, 2^2^ 0) 

-^(Ck = ^diag(5^5^ 0, {g'^+g')) , 

that these matrices are (trivially) simultaneously diagonaliz- 
able, and more importantly, their non-zero eigenvalues reside 
in distinct subspaces. This property enables us to split the 
logarithm as in (22). 

We can also verify theorem 2 for the matrices in Eqs. 
(98) and (99) by diagonalizing the gauge boson mass matrix 
■m\{tiipY^ via a rotation through the weak mixing angle 9w'- 



(g' 



{g'' 



where 



COS 

sin 



— sin 
cos Qw 

(2)\ab 



tan^vF = g I g ■ 



Therefore, the matrices m\{h)^'Y'' and m\{h^^^)ij have the 
same non-zero eigenvalues. Thus, when carrying out the sum 
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FIG. 5: A comparison of the critical temperatures as computed using 
the following methods: the standard method (solid line); the gauge- 
independent methods described in the text, derived from the full the- 
ory (dashed lines); and by performing lattice simulations (arrow). 
Note that the lattice result is much higher than the perturbative esti- 
mations and is displayed on a separate scale. 



over eigenvalues in (21) and (25), the gauge-dependent terms 
cancel. 

Finally, as outlined in section IV, the gauge-independent 
sphaleron scale v{T) is derived from the high-temperature ef- 
fective theory. In the standard model, the sphaleron scale is 
given by (79) and (30). 



A. Numerical results 

In Fig. 5, we compare critical temperatures derived follow- 
ing our gauge-independent method with the one derived using 
the conventional procedure as a function of the gauge parame- 
ter ^. In this example, we have chosen A = 0.035 correspond- 
ing to a low Higgs mass of approximately 65 GeV - close to 
the phase transition end-point according to lattice studies. 

In this plot, we display the gauge-independent results for 
Tc at two levels of approximation: at 0{H) using (102), and 
an estimate at 0{h^) with ring re-sum (60) using 



(103) 



where 



h^[V2iho,T)-^hUT) 



aVo I 



0{h^T^) . (104) 



The 0{lv'T^) subtraction represents a careful removal of that 
term so as to not double-count the contribution in Al/[^g (T). 
At present we are able to provide only a rough estimate of 
the 0(fi,2) contribution since the expression for V2{h,T) in 



140 - 




0' 

0.02 0.04 0.06 0.08 0.10 0.12 

Quartic coupling, A 

FIG. 6: Gauge-independent one-loop critical temperature as a func- 
tion of the SM Higgs quartic self-coupling. 



Ref. [45] is given in the high-T limit only. Nonetheless, we 
expect this estimate to provide an indication of magnitude of 
effects associated with higher order contributions. 

Also included are lattice results, that yield a critical tem- 
perature of 126.8 GeV, independent of ^ by construction. Our 
estimate of the higher-order contributions included in (104) 
leads to a substantially larger value of Tc, suggesting that the 
difference between the non-perturbative and 0{h) perturba- 
tive results arises in part from the omission of higher-order 
contributions. In addition, we note that the precise definition 
of Tc as obtained from the lattice studies differs from the one 
we have employed here as well as in other perturbative anal- 
yses (For a discussion of the lattice determinations, see, e.g., 
Refs. [4, 6, 7]). We speculate that part of the difference be- 
tween the lattice and perturbative results may also be due to 
this difference in definition. 

While the gauge-independent perturbative estimation of Tc 
falls below the lattice value, it is interesting that the depen- 
dence on the relevant couplings follows the trend observed 
in non-perturbative studies. To illustrate, we plot the one-loop 
Tc as a function of the Higgs quartic self-coupling A in Fig. 6. 
We observe that increasing A increases Tc in agreement with 
our qualitative expectations in (48). As we discuss shortly, this 
trend implies that the efficiency of sphaleron-induced baryon 
number washout increases with A and, thus, with the value of 
the Higgs boson mass. This trend is also observed in non- 
perturbative studies as well as in earlier gauge-dependent per- 
turbative analyses. 

We now turn our attention to the sphaleron scale, v{T), 
which we plot in Fig. 7. We observe that in the vicinity of 
the Tc obtained at 0{h) in the full theory, v{T) drops rapidly 
to zero. This behavior makes the perturbative estimate of the 
sphaleron rate at the critical temperature highly sensitive to 
small changes in Tc- Therefore, statements about the efficacy 
of baryon number preservation susceptible to large uncertain- 
ties. To illustrate, we first consider the value of this scale at 
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A=0.035 (m„«65 GeV) 
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FIG. 7: Gauge-independent sphaleron scale, v{T) in Eq. (79) as 
a function of temperature. Also shown are the gauge-independent 
0{h) (solid circle) and approximate 0{f?) (open circle) estimates 
of Tc derived from the full theory, and the lattice result (star). 

the one-loop Tc in the full theory, 95.4 GeV. We then obtain 
for the ratio 

viTc) 

4^ = 1.35, (105) 

o(h) 

implying that, at this level of approximation, sphaleron pro- 
cesses may be sufficiently quenched in the standard model to 
preserve the baryon asymmetry (see Section VI). This con- 
clusion stands in stark contrast to the conclusion one would 
reach from the results of Monte Carlo studies, which give a 
much larger Tc — well above Tq = 76.6 GeV at which v{T) 
vanishes. On the other hand, if we take the estimate of Tc at 
0{h^) using (104), we observe that it falls above To, suggest- 
ing unsuppressed sphaleron transitions that would completely 
erase the baryon asymmetry. 

In light of this strong sensitivity to higher-order contribu- 
tions in the vicinity of Tq, one should evaluate the degree to 
which a perturbative analysis is indicative. To gain insight, we 
again turn to results of Monte Carlo simulations carried out in 
Ref. [7], which found that for a range of values for A that in- 
cludes the choice used here, the sphaleron rate Fsph would be 
sufficiently quenched so as to preserve an initial baryon asym- 
metry, even though the 0{h) gauge-independent value of Tc 
falls well below the result of Monte Carlo studies. This sit- 
uation suggests that higher order contributions to Fsph may 
compensate for the larger value of Tc, resulting in stronger 
quenching of the sphaleron transitions at a given temperature 
than one might infer from perturbation theory. 

It appears, then, that one should treat with caution any con- 
clusions drawn from perturbation theory about the degree of 
baryon washout for critical temperatures in the vicinity of Tq, 
even when gauge-independence is fully maintained. We be- 
lieve, nevertheless, that perturbation theory retains some util- 
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FIG. 8: Ratio of gauge-independent sphaleron scale and critical tem- 
perature as a function of the SM Higgs quartic self-coupling. 



ity in providing guidance as to the relevant regions of param- 
eter space where baryon number preservation is likely to be 
more or less effective. To illustrate this point, we consider 
two examples. 

First, in Fig. 8 we plot the ratio v{Tc)/Tc as a function of 
A. As expected from Monte Carlo studies, this ratio decreases 
with A, implying an increase in the sphaleron rate with mu- 
The precise value of A at which baryon number preservation 
is too weak to be phenomenologically viable is uncertain - a 
point we elaborate in Section VI. Nevertheless, the trend is 
indicative and suggests that any extension of the SM scalar 
sector that reduces the effective quartic coupling at T ~ Tc is 
likely to enhance baryon number preservation. As discussed 
in Ref. [16], for example, the inclusion of new tree-level op- 
erators involving singlet degrees of freedom may lead to just 
such a weakening of Acff (Tc) while respecting the LEP lower 
bound on the SM-like Higgs boson mass. 

As a second illustration, we consider the MSSM. It has been 
thought for some time that the presence of a light right-handed 
top squark is needed in order to preserve a sufficiently large 
initial baryon asymmetry (see Ref. [9] and references therein). 
In the conventional, gauge-dependent analysis, this conclu- 
sion can be seen by considering the top squark contribution to 
the ring re-sum 



3/2' 



(106) 



where the second term cancels against the corresponding term 
in the high-T expansion of the potential. Then, the stop 
soft mass parameter m? must be negative with a magnitude 
roughly needed to cancel the stop Debye mass contribution 
Sj(r) in the first term. In this region of parameter space, the 
stop contribution to the ring potential leads to an approximate 
Tcj)'^ term in the high-T effective potential needed for a first 
order EWPT. Since the stop Yukawa coupling is 0(1), this 
contribution can be particularly pronounced. Phenomenolog- 
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FIG. 9: a. Behavior of (107) as a function of stop mass parameter 
mj. b. Broken phase free energy curve (relative to the symmetric 
phase) as a function of temperature T. A rise in the free energy 
curve (from solid to dotted line) results in a drop of Tc (indicated by 
the arrow). 



ically, the soft mass parameter to^ for the left-handed third 
generation squark doublet must be positive and relatively large 
in order to generate a phenomenologically allowed mass for 
the lightest Higgs, leaving a tachyonic right-handed stop as 
the only alternative. 

Since we have argued that the foregoing analysis is gauge- 
dependent, one may rightly question whether this scenario re- 
mains valid when a gauge-independent perturbative treatment 
is followed. We believe that it does. We reconsider (106) 
and apply the prescription (60) to construct the ring-improved 
gauge independent potential. We evaluate the resulting con- 
tribution to the difference between the energies of the elec- 
troweak symmetric and broken minima giving 



M<^min,T)-Feff(0,T) 

. T 



12n I 



(m? + y?02 + I],(r))^/^-(m? + E,(T)) 



V 3/2 



(107) 



Referring to Fig. 9a, for any > — I]j(r) the difference is 
negative. But, as m? is lowered closer to — I]j(r), the differ- 
ence in free energy rises; i.e. the free energy curve of the bro- 
ken phase is pushed up relative to that of the symmetric phase. 
As a result Tc decreases (Fig. 9b). This drop in Tc with an 
increasingly negative soft mass-squared is, in fact, observed 
in Monte Carlo studies (see, e.g.. Fig. 9 of Ref. [19]). Since 
v{Tc)/Tc correspondingly increases, the baryon asymmetry 
is more efficiently preserved. 



VI. BARYON NUMBER PRESERVATION CRITERION 

In previous sections, we have discussed at length the issue 
of gauge-dependence and its impact on the LHS of the bound 
in Eq. (6). In this section, we revisit the assumptions used to 
obtain a numeric value appearing on the RHS of the bound. In 
Appendix A, we have collected and somewhat re-organized 
previously-employed results for the BNPC in a manner that 
makes various sources of uncertainty more apparent. The ap- 
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proximate BNPC, given in (A4) is 

6 In — — > 

-InX -ln(^^) +lnZ + h\iiK. (108) 

Here, X — — In 5* characterizes the dilution of an initial 
baryon asymmetry during the transition as per (3), AiEW / tn 
gives the duration of the transition in units of the Hubble time 
tn, and k is the ratio of fluctuation determinants in (88) af- 
ter separating out the zero-mode contributions. The factor Z 
collects a number of remaining terms that enter into the per- 
turbative rate calculation: 

z = mx.. (^) , (109) 

where, A/ti and Mot are normalization factors associated with 
translational and rotational zero mode fluctuations about the 
sphaleron, Viot is the corresponding rotational volume, f is 
the number of families of fermions, and uj- is the frequency 
associated with the unstable mode of the sphaleron. Taking 
iff ft! (3 X 10-2)Mpi/r2 wifli Mpi being the Planck mass, 
and 

uj-=gvT (110) 

with being a function of T, v and the SM couplings, we 
arrive at 

4nBviTc) v{Tc) ^ 

7 m > 

g Tc Tc 

-InX -ln(^^) +lnQT+hliiK, (111) 
V iff / 

with 

Q = (3 X 10-) X (^) . 

(112) 

From Eqs. (Ill) and (112) one may attempt to derive a re- 
quirement on the ratio v{Tc)/Tc- The rotational and trans- 
lational model factors have been computed in Ref. [47] as a 
function of \/g^. These authors find that A/ji- (A/^V)i-ot ~ 
7000 and nearly constant for a substantial range in A/g^. The 
remaining terms in Eqs. (Ill) and (112) contain a stronger 
dependence on the couplings and larger theoretical uncertain- 
ties. We have explored the quantitative impact in each case 
and discuss each below. Before doing so, we adopt a repre- 
sentative set of "benchmark" values: 

T =1 , X^l 
AtEw/tH = l , Tc = 100GeV (113) 
A/.g2 = 1 , K = 10~^ 

For this choice, the inequality in Eq. (Ill) leads to 



We observe that this requirement is somewhat more relaxed 
than the one usually quoted in the literature. As we now dis- 
cuss, however, there exists considerable room for variation in 
the RHS of Eq. (114), which may be nearly a factor of two 
larger under certain circumstances. 

(i) The function T depends strongly on which fluctuations 
dominate the unstable direction, namely, those associ- 
ated with the Higgs scalar or those associated with the 
gauge degrees of freedom. The authors of Ref. [46] ar- 
gue that in the former instance one has T = 1 ["case 

(a) "] while in the latter T « {v/T^ x (l/27r) ["case 

(b) "]. Our benchmark choice corresponds to case (a). 
Generally speaking, case (b) leads to a more relaxed re- 
quirement on v{Tc)/Tc, reducing the RHS by ~ 0.1. 

(ii) The benchmark critical temperature lies between the 
one-loop, gauge-independent perturbative value and the 
value obtained from Monte Carlo studies. Increasing it 
to the lattice value lowers the RHS of Eq. (114) five 
percent or less. 

(iii) Our choice of X = 1 corresponds to allowing an ini- 
tial baryon asymmetry to be diluted by a factor of 1/e 
over the course of the EWPT. In earlier work (see, e.g., 
Ref. [9]), a value of X w 10 was used, corresponding 
to a dilution of roughly five orders of magnitude. In 
this case, the RHS of Eq. (1 14) becomes smaller by ten 
percent or more. However, recent computations of the 
initial baryon asymmetry using state-of-the art transport 
theory imply that such a dilution is likely to be wildly 
unrealistic (see, e.g., Ref. [30] and references therein). 
In the MSSM, for example, one requires a rather nar- 
row window on the relevant CP-violating parameters 
to obtain an initial baryon asymmetry that is close to 
the present value [30, 50]. A more realistic value of X 
would likely lie in the range 0.01 to 0.1, corresponding 
to a dilution of one to ten percent. Choosing X = 0.01 
increases the RHS of Eq. (1 14) by order 30%. 

(iv) According to the numerical work of Ref. [23], the value 
of the fluctuation determinant ratio k depends strongly 
on X/g^. For a Higgs mass of ^ 65 GeV, one has 
K w 10^^, which we have used as our benchmark 
value. Increasing riin to the LEP lower bound (or 
X/g^ ~ 0.3), increases k to ^ 10~^ (see Fig. 3 of 
Ref. [23]) where it is nearly maximal, while it falls back 
again to ~ lO^"' for A/g^ ^ 1. Taking k close to its 
maximal value leads to the most restrictive requirement 
on v{Tc)/Tc, raising the RHS of Eq. (114) by - 0.5. 

(v) The duration of the transition is often taken to be of or- 
der t//, and we have correspondingly used AtEw/^H ~ 
1 as our benchmark value. However, the analysis of 
Ref. [7] suggests that the transition could be as much as 
three orders of magnitude shorter In this case, the re- 
quirement on v{Tc)/Tc could be weakened by nearly 
a factor of two. 



^ > 0.73 



(114) 



In short, there exists considerable latitude in the require- 
ment for baryon number preservation. Based on the forego- 
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ing numerical exploration we take the corresponding RHS of 
Eq. (1 14) to vary between 0.4 and 1.4 as a reasonable range 
for the Standard Model. This range depends on the value of 
the scalar self-coupling(s), the duration of the transition, and 
the degree of baryon number preservation needed as dictated 
by the CP- violating transport dynamics. In any realistic analy- 
sis, one must treat all of these considerations self-consistently. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper, we have identified the origin of gauge depen- 
dence of the BNPC. We resolve these issues by (1) satisfy- 
ing degeneracy and minimization conditions at each order in 
the loop expansion and (2) identifying the gauge-independent 
sphaleron scale which arises from dimensional reduction. As 
an added bonus, our method to determine the critical temper- 
ature is substantially faster than the currently employed gauge 
dependent methods. 

After a comparison of the gauge-independent Tc in the 
SM with the results of lattice computations [51] [7], we find 
that the one-loop result appears to underestimate the value of 
the critical temperature and that a two-loop determination is 
likely to be required for a realistic perturbative approxima- 
tion. Using existing results in the literature for Kiff (0, T) at 
two-loop order in the high temperature regime, we estimate 
the numerical impact of going beyond one-loop order We em- 
phasize, however, that a robust, gauge-independent treatment 
of Tc at this order will require a future complete two-loop 
computation. When analyzing the dependence of the gauge- 
independent, perturbative computations of Tc and v{Tc) on 
the underlying parameters, we also find that the results repro- 
duce the trends obtained from Monte Carlo studies. 

Nevertheless, we believe the precise quantitative impli- 
cations of a perturbative computation for baryon number 
preservation are subject to considerable uncertainties, partic- 
ularly for Tc in the vicinity of Tq. Consequently, perturba- 
tive computations — while necessary when carrying out phe- 
nomenological investigations of the EWPT — should be taken 
as indicative rather than definitive. This need for such cau- 
tion is even more apparent when one applies the results to the 
BNPC. 

Looking to future phenomenological studies that rely on 
perturbation theory, it is apparent that previous work rely- 
ing on the use of Eq. (6) to determine the viability of EWB 
in given BSM scenarios should be revisited, taking into ac- 
count the foregoing considerations: an appropriately gauge- 
independent determination of Tc', use of a gauge-independent 
scale — analogous to v{T) — for the sphaleron rate; utilization 
of a appropriate values of X and in (111); and evaluation 
of the dependence of the fluctuation determinant k on the rel- 
evant couplings. As emphasized above, obtaining a realistic 
value of Tc is likely to require computation of the two-loop, 
finite-temperature effective potential, at least in cases such as 
the MSSM wherein the dynamics of symmetry-breaking are 
dominated by loop effects. Determining the analog of v{T) 
will necessitate obtaining the sphaleron solution to the classi- 
cal equations of motion in the presence of the extended scalar 



sector. Arriving at the appropriate value of X will rely on the 
results of quantum transport computations that yield the initial 
baryon asymmetry, while obtaining the value of F will require 
more carefully identifying the degrees of freedom that dom- 
inate the unstable mode of fluctuations about the sphaleron. 
Together with the computation of k,, this program is Ukely to 
be numerically intensive. We expect that the results will at 
least be indicative of the regions of parameter space that are 
most likely to preserve an initial baryon asymmetry, pointing 
the way to more focused and robust non-perturbative calcula- 
tions. 
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Appendix A: Approximate BNPC 

In this appendix, we summarize the steps to derive the 
washout criterion in Eqn. (6) in the standard model. We are 
generally following the discussion found in the pedagogical 
review by Quiros [9]. 

The SM expression of the baryon density ub depletion rate 
(1) is given by 



duB _ 13n/ Tsph 



dt 



2 



nB , 



(Al) 



where Uf = 3 is the number of standard model generations. 
Upon integrating, we find that the dilution at time A<ew after 
the onset of the transition at time t = is 



nB(AtEw ) 
nB(0) 



= exp 



13n 



dt 



rsph(r(t)) 

VT^{t) 



(A2) 



where T{t = 0) = T^r ~ Tc- The washout criterion corre- 
sponds to imposing a lower bound e^^ on the dilution factor 



nB{M 



ns(0) 



> e 



-X 



(A3) 



To obtain an approximate expression for the washout criterion 
we assume that the integrand in (A2) is approximately con- 
stant over the time of the transition. After taking the double 
logarithm of both sides and using (89) we find the following 
bound on v{Tc)/Tc (the baryon number preservation crite- 
rion — BNPC): 



AttB v(Tc) 
9 Tc 



6 In 



v(Tc 



Tc 
\nX ~\n 



> 



At 



EW 



tH 



\nZ + h\nK, (A4) 
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where we have chosen to normahze the duration of the phase 
transition Afgw against the Hubble time t^, and 



H 



(A5) 



contains the normahzations of the transitional and rotational 
modes of the sphaleron. As emphasized in the main text, 
the ratio v{Tc)/Tc is gauge-independent only when the 
spahleron scale v{T) arising from dimensional reduction is 
used (section IV), and when the critical temperature is com- 
puted in a manner consistent with Nielsen's identity (section 

ni). 



Appendix B: Identities in the general model 

In this appendix, we derive two useful theorems involving 
the field-dependent mass matrices in the general model, intro- 
duced in section II A. The definitions of the mass matrices are 
reproduced here for convenience: 



(Bl) 



Theorem 1. Mass matrices Mf^{(j3) and m\{<j3)ij are simul- 
tanesouly diagonalizable when fields are set equal to their 
tree-level minimum. 

Proof. Since V{^) belongs in the Lagrangian, it is invariant 
under transformations of Q, 



(B2) 



where = (1 + ct°'gT^^)ij^j is the infinitesimal transfor- 
mation law for scalar fields. Then, 

V'{<^') ^ y'((i + a'^.gr^e)*) 

= y(<i>) + a"(.gT"(l>),|^ (B3) 



We differentiate with respect to to find [31] 



dV 

gT^,-^ + {gT'^^),Mf^{^) 



0, 



(B5) 



where we have used 



a^v 



— Mfj{^). When is set equal 



to the tree-level minimum, {(l)o)i, the first term vanishes. After 
multiplying from the left by {gT"^^)k and summing over a we 
find, after using (gr"0)fc(.gr°0), = m\{q^)k^ , 

that the mass matrices are not only simultaneously diagonal- 
izable, but their eigenvalues live in distinct subspaces. □ 
Theorem 2. The gauge boson mass matrix ■m\{<j>)°''' and the 
gauge-fixing scalar boson mass matrix 'nr^{(j))ij share identi- 
cal non-zero eigenvalues. This holds for any value of(j). 

Proof. Let oj'' be an eigenvector of to^((/))°^ with eigenvalue 
A: 



Now multiply at left by {gT°-(t))j and sum over a. 

{gT-cl,),{gT-cl,UgT'c^),uj' = {gT-(P)jOj'' 

^ V ' 



(B6) 



(B7) 



= AfuJi, we have 



Then, defining {gT^( 

m\{4>)j.iNuJi = XAfujj , 



(B8) 



that is, AfuJi is an eigenvector of m\{(j))ij with the same 
eigenvalue A. 

Similarly, we could have started by letting LUi be an eigen- 
vector of m\{(f))ij with eigenvalue A. Then, by multiplying at 



left {gT''4))t, we find after defining (gT'^Cf 



the same result. 



(B9) 
□ 



implies 



dV 

(gT''$).— = 0. (B4) 



We note that although {m^Y^ and {m\)ij have the same 
non-zero eigenvalues, multiplicities of the zero eigenvalues 
are expected to be different, as they have different dimensions. 



Appendix C: High-T expansion of Kff in the general model 



In this appendix, we derive the high-tempreature expansion of the finite temperature effective potential in the general model 
and prove that the O(T^) terms (thermal masses) are gauge-independent. We start with the expression given in (25), 



scalar, i 



gauge, a 



gauge, a 



(CD 
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and use the high temperature expansion of the thermal bosonic functions given in (28) 



45 12 
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to obtain 



yeff(0,T)«Kee('/')+V^Cw(</') 
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scalars,2 

.3/2] 



scalai"s,z 



gauge, a 



scalars,2 



3 E KW]'lnK(0)/T2)- E [?m^(0)]'ln(em^(0)/T2 



(C2) 



(C3) 



where and Tig are numbers of scalar and gauge degrees of freedom as determined by the size of their respective mass matrices. 
The last three terms proportional to ln(m^/T^) combine with the Coleman-Weinberg potential Vcw(0) in (21) to cancel field- 
dependence inside the logarithms to generate terms of the form ln(r^//i^). We chose the renormalization scale /i = T so that 
these logarithmic terms vanish. The Stefan-Boltzman T'^ terms are also independent of the mean field (f). Thus, they do not 
influence the phase transition and may be dropped from the high-T expansion. Thus, at high temperature, the effective potential 
is dominated by of V;,-ee(0) and 0{T'^) and 0{T) terms. 



scalai"s,i 



Kff(0,T)«Kee(<^) + ^[ E™'(<^;0 + (3-e) E"^'(<^) 

-f E 



gauge, o 



T 
12^ 



(0;e))'/' + (3-e^/^) EK(0)) 



3/2' 



(C4) 



We demonstrate the gauge-independence of ©(T^) terms in the following way: since the sums run over the eigenvalues of the 
mass matrices, we may write them as traces of corresponding mass matrices. 

0(T2) = ^ [Tr(M2 (0) + iT^^\{^),.^ + (3 - OTrm^(0)'^' 



24 L»j 
24 



TrAf2(^)+3Trmi(</)) 



ah 



In the second line, we have used theorem 2 to cancel the gauge dependent terms, thus completing our explicit proof for gauge- 
independence of thermal masses at 0(K). 
Therefore, the high-T expansion of Kff reads 



^eff T) « Kee ('^) + ^ [ Tr (0) + 3 Tr (0)"^ 



^[ EK(<^;^))'^'+(3-^'^') EK(0)) 



3/2' 



(C5) 



scalars,? 



gauge, a 



In the language of (29), the _D-coefficient is derived from the gauge-independent 0{T'^) term; Tq is, in part, governed by the 
tree-level mass terms present in VJree('7^); the E'-coefficient comes from the non-analytic 0{T) terms; and the A-coefficient also 
comes from Vlree(^) and additional terms not displayed in (C5). 
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Appendix D: Gauge independence of tlie critical temperature 

The critical temperature is defined when two points in the 
effective potential are degenerate (26) and stationary (27). In 
this appendix we show that it follows from Nielsen's identity 
that the critical temperature is gauge-independent. 

Let Vsff{(pi,T; ^) denote the gauge-dependent finite- 
temperature effective potential. Suppose the phase transition 
takes the system from the first phase 0'^'' to the second phase 



a(2) 



at the critical temperature defined by conditions (26) 

f\T,;O=0 (Dl) 



and (27) reproduced here, 

^eff 



eff 



= 0. 



(D2) 



For now, we assume that simultaneously inverting these equa- 
tions will yield gauge dependent field-values and critical tem- 
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perature, (j)]'' = 0^*^(0, and = TdO- A total differential 
of (Dl) with respect to ^ gives 



eff 



I'eff 



(D3) 



By the stationarity condition (D2) and Nielsen's identity (33) 
the first and third terms vanish in each set of brackets, leaving 
us with 



eff 



V dT 



eff 



dT 



dT, 



0. 



(D4) 



The difference inside the brackets is, in general, non vanish- 
ing. Hence we must have 



dTc 



0, 



(D5) 



the critical temperature is gauge independent. 



dent. Start by differentiating (D2) with respect to ^. 



S^^eff 



dT, , d-'V.n 



(El) 

By (D5) the second term vanishes. We can re-express the third 
term by differentiating Nielsen's identity (33), and evaluating 
the result at (if^) and T.: 



9'^eff _ dCj dV,s ^ d^V,s 
did(f>i dcj)i d(j)j ^ d(j)idcj)j 



(E2) 



where the second term vanishes at the minimum 0^^^ . So, (El) 
becomes 



di 



(E3) 



Then, provided the curvature is non-vanishing at the 

minimum, we have 



Appendix E: Gauge dependence of scalar minimizing field ^min 

In this section, we show that the field (/)min that minimizes 
the effective potential at finite temperature is gauge depen- 



d^ 



(1) 



di 



(E4) 



that the minimizing field is inherently gauge-dependent. 
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Appendix F: One-loop effective potential in the Standard Model 

In this appendix, we provide expressions for the one-loop effective potential, with the dependence on the gauge parameter 
i shown explicitly needed for the discussion in section V. In terms of the field dependent masses in (100), the temperature- 
independent and -dependent parts of the one-loop potential are 
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where "free" represents a free-field subtraction. 
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